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Abstract 


How the initial fluctuation affects on the elliptic flow is investigated by investigating the rapidity, 


transverse 4-velocity, centrality dependencies of elliptic flow for Au+Au at 1 GeV/A with the help of an Isospin 


Quantum Molecular Dynamics (IQMD). In addition, we compare the flow calculated with respect to participant plane 


created by the initial geometry in coordinate space with the flow reconstructed by the experimental event-plane 


method, and compare the flow with the experimental data of the FOPI collaboration. It shows that there exists some 


discrepancy between the flows reconstructed by the above two methods. 
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1 Introduction 


The characterization of collective flow has proven to 
be one of the more powerful probes of the dynamics in 
heavy in collisions. Elliptic flow is an excellent 
collective flow observable which has been studied 
extensively at Bevalac and SIS as well as AGS, SPS, 
and RHIC. 
transport 


In intermediate energy, microscopic 


model calculations have stressed the 
importance of elliptic flow for extracting the equation 
of state (EOS) of nuclear matter. 

The elliptic flow is mathematically defined as 
the second coefficient of the Fourier expansion of 
particle azimuthal distribution with respect to the 
reaction plane. The origin of the elliptic flow is the 
initial geometry anisotropy of the collision system 
which originates from the uneven distribution of the 
density in the early time of the collisions. Subsequent 
dynamical evolution of the system transforms the 
anisotropy from the coordinate space to the 
momentum space, which leads to the collective motion 
observed in the final state. Based on this physical 
picture, it can be imagined that the elliptic flow should 


be very sensitive to the initial state. 


studies of the initial 


fluctuation effect on anisotropic flow in relativistic 


Recently, extensive 
nuclear collisions have been made. The initial 
fluctuation can be illustrated as the event-by-event 
fluctuation in the shape of the overlap region created 
in initial collisions. The initial eccentricity (e2), which 
is a parameter used for quantifying the initial spatial 
anisotropy, is affected by the initial spatial anisotropy. 
In high energy collisions, it has already been 
predicted by model calculations that the anisotropic 
flow considering the initial fluctuation is larger”! Han 
et al. shows that the ratio of the elliptic flow to initial 
eccentricity (v/e2) is sensitive to the initial 
fluctuation”. The value of vys) stands for the 
conversion from the initial geometry anisotropy to the 
final momentum anisotropy. Many studies prove that 
the hydrodynamic calculations on flow are in good 
agreement with the experimental data, when the initial 


[3-10] Furthermore, the 


fluctuation is taken into account 
initial fluctuation makes higher odd-order harmonic 
flow nontrivial. Previous studies have shown that the 
triangle flow is non-negligible, after the initial 
fluctuation is taken into account!!!" 

However, the initial fluctuation effect on flow 


has only been investigated in high energy collisions so 
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far. In the intermediate energy domain, its effect is not 
addressed yet. 

In this article, the initial fluctuation effect on 
the collective flow will be demonstrated in the 
intermediate energy collisions. A transport model, 
named as an Isospin Quantum Molecular Dynamics 
(IQMD), which allows the generation of events with 
event-by-event fluctuating initial conditions, is used. 
We will also compare the flow results using two 
different methods. Finally, the centrality dependence 
of the eccentricity 2 and the ratio of v2/e) will be 
presented. 


2 Analysis method 


In the collisions, the particle azimuthal distribution 


relative to the reaction plane is not isotropic and is 


usually expanded in the Fourier series'!?). 


3 2 = 
an Ea (1+9 2v, cos(n(~-Wpp))) D 
d'p 2r p,dp,dy n=l 


where the v,=<cos[n(g-wrp)|> coefficients are 
normally referred to n-th collective flow or anisotropic 
flow''*!, and the angle bracket means an average over 
all particles in all events. yrp is the reaction plane 
angle. The reaction plane is spanned by the vector of 
the impact parameter and the beam direction. It can 
not be measured experimentally, but can be estimated 


in several ways. 


2.1 Initial fluctuation effect on the reaction plane 


The original reaction plane is XZ plane in the model 
where Xpp represents for X-axis, however, the initial 
fluctuation has effects on the reaction plane as 
illustrated in Fig.1. 

The effect of the initial fluctuation on the 
reaction plane is that the initial fluctuation makes the 
participant plane (xpp) deviate from the reaction plane 
(xpp)!'*!. The participant plane angle is defined as: 
< a sin(n@) > a © 
<r“ cos(n@) > 


y, {PP} = l [arctan 
n 


where the r and ø are the coordinate position and 
azimuthal angle of each particle and <---> is the 
participant density weighted in the initial state. 

The n-th order eccentricity calculated with 
respect to this participant plane is defined as: 


<r’ cos(ng) >° + <r’ sin(ng) >? 


<r’ > (3) 


When n=2, the above two variables v and &, 


correspond to elliptic flow and eccentricity, 
respectively. In high energy heavy-ion collisions, it is 
found that the eccentricity taking the initial fluctuation 
into account is different from the eccentricity without 


the consideration of initial fluctuation!!, 


Yrp 


YPP 
h 


XRP 


Fig.1 Reaction plane and participant plane coordinate 
systems. 


2.2 Experimental method 


A common used method in experimental analysis is 


17,18 


the event-plane method!'”!8!, Tt uses the event-plane 


angle determined from the observed collective flow 
itself as an approximate reaction planel'*). 


The event plane angle is given as, 


y, {EP} = arctan (O, y O,,)/n 
Q, = > Q, COS(NG@, ) 


Q, = 2, @, sin(ng,) (4) 


in which the sum runs over all the particles used in the 
reconstruction of the event plane. The g; and œ; are the 
azimuthal angle and weight for particle i. We choose 
@i=pr for yo>0.3, @;=—pr if yo<0.3""*!. The observed 
v, with respect to the event plane is written by: 


v =< cos[n(g, - y, {EP})] > (5) 


The event plane differs in general from the 
original reaction plane for the finite multiplicity in 
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events. Therefore, the v, has to be corrected by the 
event plane resolution Rn. 

The event plane resolution for each harmonic 
is given by: 


R, =< cos[n(y,, —Wep)] > (6) 


where the angle bracket means an average over a large 
event sample. The event plane resolution depends on 
the multiplicity of particles which are used to define 
the flow vector as well as the average flow of these 


particles via the resolution parameter!!8?071) 


X=v,VM 


RO= xL INLD LID (7) 


where I is the modified Bessel function. To calculate 
the resolution we divide the full events up into two 


independent sub-events with an equal multiplicity???" 
Each sub-event resolution can be defined as: 
Rs =< coslas -yP> (8) 


where A and B denote the two subgroups of particles. 
For the given R,, the solution for y in Eq.(7) is done 
by iteration. The full event plane resolution is 
obtained by: 


Keun = RN 2 Xan) (9) 


The final collective flow measured with 
respect to the event plane can be written as: 


obs 
v 


ar (10) 


In this method, the event plane is calculated by 
the final momentum space. In this case, the flows 
reconstructed by the event plane method can be 
affected by the evolution of the dynamics. In Ref.[1] it 
shows that in Au+Au collisions at Aa =200 GeV 
from a multiphase transport model the elliptic flow v2 
and the triangle flow v3relative to the event plane are 
larger than those relative to the participant plane. 
Therefore, it is reasonable to assume that the evolution 
of the dynamics has influence on the collective flow. 


3 Results and discussion 


Figure 2 shows the yo and u dependencies of the 
elliptic flow measured with respect to the reaction 


plane ({RP}), participant plane ({PP}) and event plane 
({EP}), where yo=y/yp is the scaled rapidity and 
Uyo=Uy/Up is the scaled transverse 4-velocity. Please note 
that in the v2 calculation, we have applied the detector 
geometry cut of the FOPI detector so that a 


quantitative comparison can be made!“ 


Protons 


@ % {RP} 
a % {PP} 
+ {EP} 
* FOPI 


-0.1-0.5 © 0.5 1.0 
Yo 


Fig.2 Comparison for v,{RP} (red circle), v,{PP} (blue 
triangle) and v2{EP} (green diamond) in Au+Au collisions at 
E=1 GeV/A and centrality 0.25<b<0.45 by the IQMD model 
with a parameter set of Soft and Momentum-dependent 
equation of state (SM). The black stars are the experimental 
data from the FOPI Collaboration. The panel (a) plots yo 
dependence where v, is integrated over uw, but constrained to 
uy>0.4. The panel (b) shows the uo dependence in the indicated 
yo bin. 

The rapidity dependence of Fig.2(a) shows a 
V-shape. The proton shows an in-plane emission (v2>0) 
in project-like and target-like regions (larger rapidity 
value in absolute value) while it shows a squeeze-out 
emission in the overlapping region (mid-rapidity). 
From the mid-rapidity to the projectile-like/target-like 
rapidities, the proton favors the squeeze-out emission 
to the in-plane emission. This is consistent with the 
shadowing effect. 

The transverse velocity dependence (Fig.2(b)) 
shows that with the increase of wo, the squeeze-out 
emission is stronger, i.e. the proton with higher 
velocity can be more easily escaped from overlapping 
zone. 

From Fig.2, we can see that the value of v.{RP} 
is a little larger than v.{PP}. However, there is almost 
no discrepancy between v2{EP} and v2{PP}. That 
means vz becomes weak by the initial fluctuation but is 
not sensitive to the dynamical evolution. This 
phenomenon is different from what is known at high 
energies, where v2 is both enlarged by the initial 
fluctuation and the dynamical evolution. In the figure, 
the FOPI data is also plotted in order to check our 
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model calculations. Even though quantitative 


agreements of various calculations are not reached, the 
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trend of vz as functions of yo and uwis consistent with 
the data. 
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Fig.3 v of protons in Au+Au collisions for different centrality ranges. In the upper three panels, v, is calculated for transverse 
4-velocities u,>0.4. In the lower three panels, v2 is calculated for the rapidity |yo|<0.4. 
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Fig. 4 & and v»/é as a function of impact parameter for 
Au+Au collisions at E=1GeV/A. 

Figure 3 displays that vz versus yo and uw in 
three different centralities. While the shapes look 
similar in all centralities, the elliptic flow of protons in 
mid-rapidity (upper panels) or high velocity (lower 
panels) regions shows stronger squeeze-out effect at 


some intermediate centralities (0.45<bo<0.55). With 
the increase of bo, the initial fluctuation and dynamical 
evolution also play an increasing role as seen from the 
behavior of increasing difference between v2{PP}, 
v2{RP} and v2{EP}. Again, the similarity between 
v2{PP} and v.{EP} indicates that there is almost no 
effect by the evolution of dynamics. 

The upper panel of Fig.4 shows the impact 
parameter (b) dependence of eccentricity e2, It shows 
that £2 increases with b and the ¢2{PP} is larger than 
é.{EP}. The lower panel shows the ratio v2/ée. as a 
function of b. We can see the absolute values of the 
ratio also increase with b and |v. {PP}/e, | is smaller 
than |v.{RP}/g|. The smaller the absolute ratio is, the 
less efficient the conversion from the initial geometry 
anisotropy to the final momentum anisotropy is. 


4 Conclusion 


The elliptic flow of Aut+Au collision at 1 GeV/A by 
three different methods for the determination of 
reaction plane was studied. And it was found that the 
initial event-by-event geometry fluctuation had effects 
on the final-state elliptic flow. In the present 
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calculation, the initial fluctuation weakens the elliptic 
flow, in other words the initial fluctuation makes the 
conversion from initial spatial eccentricity to final 
momentum smaller. In this paper, we studied the 
dependence of é¢ on impact parameter with the 
consideration of the initial fluctuation and compared 
the ratio v2/e. with or without the consideration of the 
initial fluctuation. It shows that the initial fluctuation 
enhances the eccentricity e2 but it decreases the ratio 
of |v2/é2|. Moreover, we also found that the dynamical 
evolution has little influence on the elliptic flow in our 
studied cases. Our simulation gives a similar trend for 
rapidity and velocity dependencies, but the results are 
quantitatively smaller than the experimental data for 
mid-rapidity and high velocity protons, which may be 
caused by two reasons. One is the variation of physical 
parameters (such as the ground state densities, 
interaction ranges) whose precise values are not 
known; the other is the consequence of the 
impossibility to build a ground state nucleus with all 
its detailed structure in a semi-classical molecular 


dynamics approach”), 
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